Efficient identification of QTL for agronomic traits in foxtail millet (Setaria italica) using RTM- and MLM-GWAS.

KEY MESSAGE: Eleven QTLs for agronomic traits were identified by RTM- and MLM-GWAS, putative candidate genes were predicted and two markers for grain weight were developed and validated. Foxtail millet (Setaria italica), the second most cultivated millet crop after pearl millet, is an important grain crop in arid regions. Seven agronomic traits of 408 diverse foxtail millet accessions from 15 provinces in China were evaluated in three environments. They were clustered into two divergent groups based on genotypic data using ADMIXTURE, which was highly consistent with their geographical distribution. Two models for genome-wide association studies (GWAS), namely restricted two-stage multi-locus multi-allele (RTM)-GWAS and mixed linear model (MLM)-GWAS, were used to dissect the genetic architecture of the agronomic traits based on 13,723 SNPs. Eleven quantitative trait loci (QTLs) for seven traits were identified using two models (RTM- and MLM-GWAS). Among them, five were considered stable QTLs that were identified in at least two environments using MLM-GWAS. One putative candidate gene (SETIT_006045mg, Chr4: 744,701-746,852) that can enhance grain weight per panicle was identified based on homologous gene comparison and gene expression analysis and was validated by haplotype analysis of 330 accessions with high-depth (10×) resequencing data (unpublished). In addition, homologous gene comparison and haplotype analysis identified one putative foxtail millet ortholog (SETIT_032906mg, Chr2: 5,020,600-5,029,771) with rice affecting the target traits. Two markers (cGWP6045 and kTGW2906) were developed and validated and can be used for marker-assisted selection of foxtail millet with high grain weight. The results provide a fundamental resource for foxtail millet genetic research and breeding and demonstrate the power of integrating RTM- and MLM-GWAS approaches as a complementary strategy for investigating complex traits in foxtail millet.

Combined linkage analysis and association mapping identifies genomic regions associated with yield-related and drought-tolerance traits in wheat (Triticum aestivum L.).

KEY MESSAGE: Combined linkage analysis and association mapping identified genomic regions associated with yield and drought tolerance, providing information to assist breeding for high yield and drought tolerance in wheat. Wheat (Triticum aestivum L.) is one of the most widely grown food crops and provides adequate amounts of protein to support human health. Drought stress is the most important abiotic stress constraining yield during the flowering and grain development periods. Precise targeting of genomic regions underlying yield- and drought tolerance-responsive traits would assist in breeding programs. In this study, two water treatments (well-watered, WW, and rain-fed water stress, WS) were applied, and five yield-related agronomic traits (plant height, PH; spike length, SL; spikelet number per spike, SNPS; kernel number per spike, KNPS; thousand kernel weight, TKW) and drought response values (DRVs) were used to characterize the drought sensitivity of each accession. Association mapping was performed on an association panel of 304 accessions, and linkage analysis was applied to a doubled haploid (DH) population of 152 lines. Eleven co-localized genomic regions associated with yield traits and DRV were identified in both populations. Many previously cloned key genes were located in these regions. In particular, a TKW-associated region on chromosome 2D was identified using both association mapping and linkage analysis and a key candidate gene, TraesCS2D02G142500, was detected based on gene annotation and differences in expression levels. Exonic SNPs were analyzed by sequencing the full length of TraesCS2D02G142500 in the association panel, and a rare haplotype, Hap-2, which reduced TKW to a lesser extent than Hap-1 under drought stress, and the Hap-2 varieties presented drought-insensitive. Altogether, this study provides fundamental insights into molecular targets for high yield and drought tolerance in wheat.

SiMYBS3, Encoding a Setaria italica Heterosis-Related MYB Transcription Factor, Confers Drought Tolerance in Arabidopsis.

Drought is a major limiting factor affecting grain production. Drought-tolerant crop varieties are required to ensure future grain production. Here, 5597 DEGs were identified using transcriptome data before and after drought stress in foxtail millet (Setaria italica) hybrid Zhangza 19 and its parents. A total of 607 drought-tolerant genes were screened through WGCNA, and 286 heterotic genes were screened according to the expression level. Among them, 18 genes overlapped. One gene, Seita.9G321800, encoded MYBS3 transcription factor and showed upregulated expression after drought stress. It is highly homologous with MYBS3 in maize, rice, and sorghum and was named SiMYBS3. Subcellular localization analysis showed that the SiMYBS3 protein was located in the nucleus and cytoplasm, and transactivation assay showed SiMYBS3 had transcriptional activation activity in yeast cells. Overexpression of SiMYBS3 in Arabidopsis thaliana conferred drought tolerance, insensitivity to ABA, and earlier flowering. Our results demonstrate that SiMYBS3 is a drought-related heterotic gene and it can be used for enhancing drought resistance in agricultural crop breeding.

Resistance of QYm.nau-2D to wheat yellow mosaic virus was derived from an alien introgression into common wheat.

KEY MESSAGE: The QYm.nau-2D locus conferring wheat yellow mosaic virus resistance is an exotic introgression and we developed 11 diagnostic markers tightly linked to QYm.nau-2D. Wheat yellow mosaic virus (WYMV) is a serious disease of winter wheat in China. Breeding resistant varieties is the most effective strategy for WYMV control. A WYMV resistant locus QYm.nau-2D on the chromosome arm 2DL has been repeatedly reported but the mapped region is large. In the present study, we screened recombinants using a biparental population and mapped QYm.nau-2D into an 18.8 Mb physical interval. By genome-wide association studies of 372 wheat varieties for WYMV resistance in four environments, we narrowed down QYm.nau-2D into a 16.4 Mb interval. Haplotype analysis indicated QYm.nau-2D were present as six different states due to recombination during hybridization breeding. QYm.nau-2D was finally mapped into a linkage block of 11.2 Mb. Chromosome painting using 2D specific probes and collinearity analysis among the published sequences corresponding to QYm.nau-2D region indicated the block was an exotic introgression. The Illumina-sequenced reads of four diploid Aegilops species were mapped to the sequence of Fielder, a variety having the introgression. The mapping reads were significantly increased at the putative introgression regions of Fielder. Ae. uniaristata (NN) had the highest mapping reads, suggesting that QYm.nau-2D was possibly an introgression from genome N. We investigated the agronomic performances of different haplotypes and observed no linkage drag of the alien introgression for the 15 tested traits. For marker-assisted selection of QYm.nau-2D, we developed 11 diagnostic markers tightly linked to the locus. This research provided a case study of an exotic introgression, which has been utilized in wheat improvement for WYMV resistance.

Identification of Brachypodium distachyon B3 genes reveals that BdB3-54 regulates primary root growth.

B3 is a class of plant-specific transcription factors with important roles in plant development and growth. Here, we identified 69 B3 transcription factors in Brachypodium distachyon that were unevenly distributed across all five chromosomes. The ARF, REM, LAV, and RAV subfamilies were grouped based on sequence characteristics and phylogenetic relationships. The phylogenetically related members in the B3 family shared conserved domains and gene structures. Expression profiles showed that B3 genes were widely expressed in different tissues and varied in response to different abiotic stresses. BdB3-54 protein from the REM subfamily was located in the nucleus by subcellular localization and processed transcriptional activation activity. Overexpression of BdB3-54 in Arabidopsis increased primary root length. Our study provides a basis for further research on the functions of BdB3 genes.

Identification of GATA Transcription Factors in Brachypodium distachyon and Functional Characterization of BdGATA13 in Drought Tolerance and Response to Gibberellins.

GATA transcription factors (TFs) are type IV zinc-finger proteins that have roles in plant development and growth. The 27 GATA TFs identified in the Brachypodium distachyon genome in this study were unevenly distributed across all five chromosomes and classified into four subgroups. Phylogenesis-related GATAs shared similar gene structures and conserved motifs. Expression profiles showed that all BdGATA genes were expressed in leaves and most were induced by PEG treatment. BdGATA13 was predominantly expressed in leaf tissue and phylogenetically close to OsSNFL1, AtGNC, and AtGNL. Its protein was detected in the nucleus by subcellular localization analysis. Overexpression of BdGATA13 in transgenic Arabidopsis resulted in darker green leaves, later flowering, and more importantly, enhanced drought tolerance compared to the wild type. BdGATA13 also promoted primary root development under GA treatment. These results lay a foundation for better understanding the function of GATA genes in B. distachyon and other plants.

Dissection and cytological mapping of chromosome arm 4VS by the development of wheat-Haynaldia villosa structural aberration library.

KEY MESSAGE: A cytological map of Haynaldia villosa chromosome arm 4VS was constructed to facilitate the identification and utilization of beneficial genes on 4VS. Induction of wheat-alien chromosomal structure aberrations not only provides new germplasm for wheat improvement, but also allows assignment of favorable genes to define physical regions. Especially, the translocation or introgression lines carrying alien chromosomal fragments with different sizes are useful for breeding and alien gene mapping. Chromosome arm 4VS of Haynaldia villosa (L.) Schur (syn. Dasypyrum villosum (L.) P. Candargy) confers resistances to eyespot and wheat yellow mosaic virus (WYMV). In this research, we used both irradiation and the pairing homoeologous gene (Ph) mutant to induce chromosomal aberrations or translocations. By using the two approaches, a structural aberration library of chromosome arm 4VS was constructed. In this library, there are 57 homozygous structural aberrations, in which, 39 were induced by the Triticum aestivum cv. Chinese Spring (CS) ph1b mutant (CS ph1b) and 18 were induced by irradiation. The aberrations included four types, i.e., terminal translocation, interstitial translocation, deletion and complex structural aberration. The 4VS cytological map was constructed by amplification in the developed homozygous aberrations using 199 4VS-specific markers, which could be allocated into 39 bins on 4VS. These bins were further assigned to their corresponding physical regions of chromosome arm 4DS based on BLASTn search of the marker sequences against the reference sequence of Aegilops tauschii Cosson. The developed genetic stocks and cytological map provide genetic stocks for wheat breeding as well as alien gene tagging.

Development and characterization of a complete set of Triticum aestivum-Roegneria ciliaris disomic addition lines.

KEY MESSAGE: A complete set wheat-R. ciliaris disomic addition lines (DALs) were characterized and the homoeologous groups and genome affinities of R. ciliaris chromosomes were determined. Wild relatives are rich gene resources for cultivated wheat. The development of alien addition chromosome lines not only greatly broadens the genetic diversity, but also provides genetic stocks for comparative genomics studies. Roegneria ciliaris (genome ScScYcYc), a tetraploid wild relative of wheat, is tolerant or resistant to many abiotic and biotic stresses. To develop a complete set of wheat-R. ciliaris disomic addition lines (DALs), we undertook a euplasmic backcrossing program to overcome allocytoplasmic effects and preferential chromosome transmission. To improve the efficiency of identifying chromosomes from Sc and Yc, we established techniques including sequential genomic in situ hybridization/fluorescence in situ hybridization (FISH) and molecular marker analysis. Fourteen DALs of wheat, each containing one pair of R. ciliaris chromosomes pairs, were characterized by FISH using four repetitive sequences [pTa794, pTa71, RcAfa and (GAA)10] as probes. One hundred and sixty-two R. ciliaris-specific markers were developed. FISH and marker analysis enabled us to assign the homoeologous groups and genome affinities of R. ciliaris chromosomes. FHB resistance evaluation in successive five growth seasons showed that the amphiploid, DA2Yc, DA5Yc and DA6Sc had improved FHB resistance, indicating their potential value in wheat improvement. The 14 DALs are likely new gene resources and will be phenotyped for more agronomic performances traits.

Sequencing flow-sorted short arm of Haynaldia villosa chromosome 4V provides insights into its molecular structure and virtual gene order.

BACKGROUND: Haynaldia villosa (H. villosa) has been recognized as a species potentially useful for wheat improvement. The availability of its genomic sequences will boost its research and application. RESULTS: In this work, the short arm of H. villosa chromosome 4V (4VS) was sorted by flow cytometry and sequenced using Illumina platform. About 170.6 Mb assembled sequences were obtained. Further analysis showed that repetitive elements accounted for about 64.6% of 4VS, while the coding fraction, which is corresponding to 1977 annotated genes, represented 1.5% of the arm. The syntenic regions of the 4VS were searched and identified on wheat group 4 chromosomes 4AL, 4BS, 4DS, Brachypodium chromosomes 1 and 4, rice chromosomes 3 and 11, and sorghum chromosomes 1, 5 and 8. Based on genome-zipper analysis, a virtual gene order comprising 735 gene loci on 4VS genome was built by referring to the Brachypodium genome, which was relatively consistent with the scaffold order determined for Ae. tauschii chromosome 4D. The homologous alleles of several cloned genes on wheat group 4 chromosomes including Rht-1 gene were identified. CONCLUSIONS: The sequences provided valuable information for mapping and positional-cloning genes located on 4VS, such as the wheat yellow mosaic virus resistance gene Wss1. The work on 4VS provided detailed insights into the genome of H. villosa, and may also serve as a model for sequencing the remaining parts of H. villosa genome.

Development of intron targeting (IT) markers specific for chromosome arm 4VS of Haynaldia villosa by chromosome sorting and next-generation sequencing.

BACKGROUND: Haynaldia villosa (L.) Schur (syn. Dasypyrum villosum L. Candargy, 2n = 14, genome VV) is the tertiary gene pool of wheat, and thus a potential resource of genes for wheat improvement. Among other, wheat yellow mosaic (WYM) resistance gene Wss1 and a take-all resistance gene were identified on the short arm of chromosome 4 V (4VS) of H. villosa. We had obtained introgressions on 4VS chromosome arm, with the objective of utilizing the target genes. However, monitoring these introgressions has been a daunting task because of inadequate knowledge as to H.villosa genome, as reflected by the lack of specific markers. RESULTS: This study aims to develop 4VS-specific markers by combination of chromosome sorting and next-generation sequencing. The short arm of chromosome 4VS of H.villosa was flow-sorted using a FACSVantage SE flow cytometer and sorter, and then sequenced by Illumina sequencing. The sequence of H. villosa 4VS was assembled by the software Hecate, and then was compared with the sequence assemblies of wheat chromosome arms 4AL, 4BS and 4DS and Ae. tauschii 4DS, with the objectives of identifying exon-exon junctions and localizing introns on chromosome 4VS of H. villosa. The intron length polymorphisms suitable for designing H. villosa primers were evaluated with criteria. Consequently, we designed a total of 359 intron targeting (IT) markers, among which 232 (64.62%) markers were specific for tracing the 4VS chromatin in the wheat background. CONCLUSION: The combination of chromosome sorting and next-generation sequencing to develop specific IT markers for 4VS of H. villosa has high success rate and specificity, thus being applicable for the development of chromosome-specific markers for alien chromatin in wheat breeding.